Radio resource management in metropolitan area networks

ABSTRACT

A system for allocating resources in a communications network includes a plurality of access points, where each access point in the plurality of access points is in communication with one or more user devices. The system also includes a central controller in communication with the plurality of access points. The central controller is configured to receive traffic information and channel information for access point—user device links in the communications network. The central controller is also configured to receive an intensity of independent homogeneous Poisson traffic intended for the one or more user devices. The central controller is further configured to determine a joint resource allocation based at least in part on the received traffic and channel information and the intensity of independent homogeneous Poisson traffic, where the joint resource allocation is determined with an algorithm.

CROSS-REFERENCE TO RELATED APPLICATION

The present application is a Continuation Application of U.S. patent application Ser. No. 16/146,276 filed on Sep. 28, 2018, which claims the priority benefit of U.S. Provisional Patent Application No. 62/565,681 filed on Sep. 29, 2017, the entire disclosures of which are incorporated herein by reference.

BACKGROUND

Recent years have witnessed an explosion of mobile data traffic due to proliferation of smart terminals (e.g., smart phones and tablet personal computers) and increasing use of mobile applications. It is predicted that global mobile data traffic will increase eightfold between 2015 and 2020. Wireless operators are investing heavily in infrastructure, including increasingly denser deployment of access points (APs) to improve cellular coverage and capacity for homes and businesses in urban, suburban, as well as rural areas.

SUMMARY

An illustrative system for allocating resources in a communications network includes a plurality of access points and a central controller in communication with the plurality of access points. Each access point in the plurality of access points is in communication with one or more user devices. The central controller is configured to receive traffic information and channel information for access point—user device links in the communications network. The central controller is also configured to determine a joint resource allocation based at least in part on the received traffic and channel information and at least in part on an algorithm. The central controller is further configured to transmit the joint resource allocation to the plurality of access points, where the joint resource allocation comprises spectrum allocation, user association, and power allocation.

An illustrative method for allocating resources in a communications network includes receiving, by a central controller and from a plurality of access points in communication with the central controller, traffic information and channel information for access point—user device links in a communications network. Each access point in the plurality of access points is in communication with one or more user devices. The method also includes determining a joint resource allocation based at least in part on the received traffic and channel information and at least in part on an algorithm. The method further includes transmitting the joint resource allocation to the plurality of access points, where the joint resource allocation comprises spectrum allocation, user association, and power allocation.

Other principal features and advantages of the invention will become apparent to those skilled in the art upon review of the following drawings, the detailed description, and the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

Illustrative embodiments of the invention will hereafter be described with reference to the accompanying drawings, wherein like numerals denote like elements.

FIG. 1 depicts a network in accordance with an illustrative embodiment.

FIG. 2 depicts a table with parameter configurations used to test the performance gain of the proposed allocation scheme in accordance with an illustrative embodiment.

FIG. 3 depicts a test deployment in a medium scale network in accordance with an illustrative embodiment.

FIG. 4 depicts the objective of average packet delay versus traffic intensity curves for the deployment of FIG. 3 in accordance with an illustrative embodiment.

FIG. 5 depicts a test deployment in a large scale network in accordance with an illustrative embodiment.

FIG. 6 depicts the objective of average packet delay versus traffic intensity curves for the deployment of FIG. 5 in accordance with an illustrative embodiment.

FIG. 7 is a block diagram depicting a central controller for the proposed system in accordance with an illustrative embodiment.

DETAILED DESCRIPTION

Described herein is a centralized radio resource management system for metropolitan area networks with a large number of access points (APs) and user equipments (UEs) or user devices. Conceived for a moderate timescale (seconds to minutes), a non-convex infinite-size problem is first formulated to allow for highly flexible spectrum allocation, user association, and power control according to large-scale channel and traffic conditions. A scalable reformulation of the problem is obtained by exploiting the sparse structure of the optimal solution. An efficient algorithm is also proposed to solve the reformulated problem with guranteed convergence. Moreover, each iteration is performed in closed form, which makes centralized resource allocation practical for networks with thousands of APs and UEs. Numerical results show that the proposed scheme yields a significantly improved throughput region and a reduced average packet delay as compared to conventional schemes.

Wireless systems have emerged as a ubiquitous part of modern data communication networks. Demand for these systems continues to grow as applications involving both voice and data expand beyond their traditional wireline service requirements. In order to meet the increasing demand in data rates, a large number of access points (APs) have been deployed to provide improved cellular coverage, capacity, and applications for homes and businesses, as well as metropolitan and rural public spaces. These APs may be densely deployed and may also utilize additional spectrum resources, such as unlicensed spectrum resources, high-frequency spectrum resources, etc. Each AP can be in the form of a computing device that includes a processor, a memory, a transceiver, etc.

Due to limited resources in current wireless networks, efficient resource allocation (e.g., spectrum allocation, power control, link scheduling, routing, and congestion control) is important in achieving high performance and providing satisfactory quality-of-service (QoS). Conventional spectrum allocation schemes include full spectrum reuse, partial frequency reuse, and fractional frequency reuse. However, the spectrum is either under-utilized or over-utilized under these traditional schemes, resulting in either low spectral efficiency or strong link interference. Moreover, these schemes are not traffic-aware in that they are not adaptive to different traffic conditions. Considering the user association schemes used in the industry, most of the existing schemes are based on simple heuristics which waste precious system capacity and lead to suboptimal performance, while offering no performance guarantees. For example, by default, in today's cellular and Wi-Fi networks users simply associate with the base station from which they receive the strongest signal. Power control is another way of mitigating inter-cell interference as well as saving energy. Instead of letting each AP always transmit using full power, some enhanced inter-cell interference coordination (eICIC) techniques, such as the almost blank subframe (ABS) control, have been proposed in long term evolution (LTE) and long term evolution advanced (LTE-A) networks. However, these types of distributed algorithms often compromise global optimality for the overall network.

Resource allocation in wireless networks has been extensively studied over the last few decades. Spectrum allocation schemes are basically designed either from an economic perspective or from a technical perspective As for the user association problem, NP-hard non-convex integer programming problems are often formulated. However, these problems make it difficult to obtain a local optimum, and extremely difficult to obtain a global optimum. Similar to the user association problem, optimal power control can also be a hard non-convex optimization problem due to the complicated interference coupling between links. Most power control algorithms are designed in a distributed manner. This is especially the case for the power control problem in a large-scale network with hundreds or even thousands of APs. However, this is at the sacrifice of system performance since distributed algorithms typically do not achieve the global optimum.

The present system addresses the preceding problems by proposing a centralized optimization-based framework for resource allocation. Unlike traditional systems, the proposed framework is able to jointly solve the problems of spectrum allocation, user association, and power control for networks with thousands of APs. An efficient algorithm with low complexity is proposed herein to obtain a near-optimal allocation solution. One key distinction of the problem formulation described herein is that, instead of assuming all APs always transmit using full power, each AP is allowed to apply a continuous power spectral density over the entire band, i.e., the size of the problem of power allocation at each AP has already grown to infinity.

The present system is the first to propose centralized joint spectrum allocation, user association, and power allocation for very large networks with thousands of APs and UEs. This is accomplished in part by formulation of a non-convex infinite-size joint spectrum allocation, user association, and power allocation problem (considering any arbitrary concave utility function), and reformulating this problem into an optimization problem with finite size which is amenable to a highly efficient algorithm. In order to solve the reformulated problem, an efficient algorithm which is guranteed to converge is proposed. Moreover, each iteration is performed in closed form, which makes centralized resource allocation practically feasible, even for a very large network. As discussed in more detail below, numerical results show that the proposed resource allocation scheme yields a significantly improved throughput region and reduced average packet delay as compared with conventional schemes.

In an illustrative embodiment, a network can include n APs and k UEs. FIG. 1 depicts a network 100 in accordance with an illustrative embodiment. As depicted, the network includes n APs and k UEs, where n and k can be any number. In an illustrative embodiment, each of the APs is a computing device that includes a processor, a memory, a transceiver, an interface, etc. The APs can also include software/algorithms stored on the memory and used to perform any of the AP operations described herein. Each of the APs is in communication with a plurality of UEs. Each of the UEs can be a smart phone, tablet, or other network device that includes a processor, memory, transceiver, interface, etc.

The network 100 also a central controller 105 in communication with the APs. The central controller 105 can also be a computing device that includes a processor, memory, transceiver, and interface. Software and/or algorithms can be stored on the memory of the central controller 105 and used to implement any of the central controller operations described herein. In an illustrative embodiment, the central controller 105 is in communication with all APs, and collects traffic and channel/interference for all AP-UE links over time. The central controller 105 can also be informed of the intensity of independent homogeneous Poisson traffic intended for every UE by the APs. Given a spectrum resource of a certain bandwidth, one task of the central controller 105 is to allocate resources in order to maximize the long-term network utility. The allocated resources include, but are not limited to, power allocation, spectrum allocation, and user association.

In one embodiment, resource allocation can be performed on a slow timescale, which is conceived to be once every few seconds or minutes based on current computational power available to most devices. This makes information exchange (channel state information feedback) and joint resource allocation viable at the central controller 105. In addition, since a period of slow-timescale resource allocation is much longer than the channel coherence time, the average channel conditions can be accurately modeled and measured using path loss and the statistics of small scale fading. The frequency resources are assumed to be homogeneous on a slow timescale.

The set of AP indexes can be denoted by N={1, . . . , n} and the set of UE indexes denoted by K={1, . . . , k}. Arbitrary association is allowed such that each AP can simultaneously serve any subset of UEs, and each UE can be simultaneously served by any subset of APs. The power constraint can be on the (transmit) power spectral density (PSD) of individual APs. The spectrum allocation is integrated into the frequency-selective power allocation, as choosing a subset of APs to transmit over a certain part of the spectrum is equivalent to allocating positive transmit power to these APs and zero power to others over that part of the spectrum.

When a receiver decodes information from a transmitter, it can be assumed that the receiver treats signals from all other transmitters as additive white Gaussian noise. Allocating more spectrum and/or power to a given AP-UE link will thus improve the throughput of this given link at the expense of additional interference to other links. The allocation problem then can be viewed as determining the PSDs of all AP-UE links over the entire spectrum in order to maximize the network utility. For convenience, the entire frequency band can be thought of as one unit. The PSD of a link from AP i to UE j can be denoted as p_(i→j)(f), f ∈[0,1], and the PSDs of all links can be denoted by a vector p(f)=(p_(i→j)(f))_(i∈N,j∈K). It can also be assumed that an AP serves at most one UE over each slice of the spectrum, resulting in the following constraint:

$\begin{matrix} \begin{matrix} {{{\sum\limits_{j \in K}{{p_{i\rightarrow j}(f)}}_{0}} \leq 1},} & {{\forall{i \in N}},{\forall{f \in {\left\lbrack {0,1} \right\rbrack.}}}} \end{matrix} & {{Eq}.\mspace{14mu} 1} \end{matrix}$

Nevertheless, each AP can serve multiple UEs over different slices of the spectrum. Assuming Eq. 1 is satisfied, power constraints on the APs can be written as:

$\begin{matrix} \begin{matrix} {{{p_{i\rightarrow j}(f)} \geq 0},} & {{\forall{i \in N}},{j \in K},{\forall{f \in {\left\lbrack {0,1} \right\rbrack.}}}} \end{matrix} & {{Eq}.\mspace{14mu} 2} \\ {\begin{matrix} {{\left( {\sum\limits_{j \in K}{\int_{0}^{1}{\left\lbrack {p_{i\rightarrow j}(f)} \right\rbrack^{\alpha}f}}} \right)^{\frac{1}{\alpha}} \leq P_{\max}},} & {{\forall{i \in N}},{\forall{f \in \left\lbrack {0,1} \right\rbrack}}} \end{matrix}.} & {{Eq}.\mspace{14mu} 3} \end{matrix}$

Equation 3 can be generalized to many different types of power constraints by tuning different α≥1. For example, when α=1, it constrains the total power of each AP over the entire spectrum, and when α→∞, it constrains the PSD of each AP over each hertz. Given p(f), the service rate to UE j contributed by AP i can be calculated as:

r _(i→j)=∫₀ ⁻¹ s(p(f))f,   Eq. 4:

where s( ) is the general instantaneous rate function.

For purposes of the network system described herein, s( )can be written as a function involving fractional (or ratio) terms, as indicated in Eq. 5 below:

$\begin{matrix} {{{s\left( {p(f)} \right)} = {q\left( \frac{A\left( {p(f)} \right)}{B\left( {p(f)} \right)} \right)}},} & {{Eq}.\mspace{14mu} 5} \end{matrix}$

where A( ) is a nonnegative function and B() is a positive function of p(f). This is generally the case because the instantaneous rate depends on the signal-to-interference-plus-noise ratio. The total service rate of UE j, which is denoted as r_(j), is calculated by summing over all APs over the entire spectrum as follows:

$\begin{matrix} {{r_{j} = {\sum\limits_{i \in N}r_{i\rightarrow j}}}.} & {{Eq}.\mspace{14mu} 6} \end{matrix}$

The fundamental resource allocation problem is to maximize the long-term utility by adapting the PSDs of all AP-UE links. Collecting the constraints of Eq. 1-6 above, P0 can be formulated as follows:

$\begin{matrix} \begin{matrix} {2\underset{r,{\{{p{(f)}}\}}}{maximize}} & {U(r)} \end{matrix} & {P0a} \\ {{{{subjectto}\mspace{14mu} r_{j}} = {\sum\limits_{i \in N}{\int_{0}^{1}{{q\left( \frac{A\left( {p(f)} \right)}{B\left( {p(f)} \right)} \right)}f}}}},{\forall{j \in K}}} & {P0b} \\ {{{\sum\limits_{j \in K}{{p_{i\rightarrow j}(f)}}_{0}} \leq 1},{\forall{i \in N}},{\forall{f \in \left\lbrack {0,1} \right\rbrack}}} & {P0c} \\ {{\left( {\sum\limits_{j \in K}{\int_{0}^{1}{\left\lbrack {p_{i\rightarrow j}(f)} \right\rbrack^{\alpha}f}}} \right)^{\frac{1}{\alpha}} \leq P_{\max}},{\forall{i \in N}},{\forall{f \in {\left\lbrack {0,1} \right\rbrack.}}}} & {P0d} \\ {{{p_{i\rightarrow j}(f)} \geq 0},{\forall{i \in N}},{j \in K},{\forall{f \in \left\lbrack {0,1} \right\rbrack}},} & {P0e} \end{matrix}$

where U(r) is the network utility function. It is assumed herein that U( ) is a concave function. In fact, most conventional utility functions are concave such as the sum rate, the minimum user service rate (max-min fairness), and the sum log-rate (proportional fairness).

It is noted that P0 can be efficiently solved using the proposed methods described herein. The methods for solving P0 are elaborated upon below by narrowing down to specific functions for concreteness. Average (negative) packet delay is used as the network utility function:

$\begin{matrix} {{{u(r)} = {- {\sum\limits_{j \in K}\frac{\lambda_{j}}{\left( {r_{j} - \lambda_{j}} \right)^{+}}}}},} & {{Eq}.\mspace{14mu} 7} \end{matrix}$

where is the homogeneous Poisson packet arrival rate of UE j, and the extended real-valued function 1/x⁻=1/x if x>0 and 1/x⁺=+∞ if x≤0 . It can be seen that 1/x⁺ is convex on (−∞, +∞). The choice of this utility function also assumes exponential packet length and a relatively conservative rate. If r_(j)≤λ, the packet delay is infinite, and the system becomes unstable.

For the instantaneous rate function, the Shannon formula can be used in one embodiment. That is

$\begin{matrix} {{r_{i\rightarrow j} = {\int_{0}^{1}{{\log\left( {1 + \frac{{p_{i\rightarrow j}(f)}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}{p_{l\rightarrow q}(f)}}}}}} \right)}f}}},} & {{Eq}.\mspace{14mu} 8} \end{matrix}$

where g_(i>1) is the average gain of link i→j which captures the effects of path loss and shadowing, and

$\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}{p_{l\rightarrow q}(f)}}}$

is the interference received from other APs operating over the same frequency f.

For the power constraints, the instantaneous power constraint is chosen by letting α→∞ in Eq. 3, as follows:

0≤p_(i→j)(f)≤P_(max), ∀i ∈N, j ∈K, ∀f ∈[0,1].   Eq. 9:

P0 can then be rewritten as P1, as follows:

$\begin{matrix} {2\underset{r,{\{{p{(f)}}\}}}{maximize}\mspace{14mu} {U(r)}} & {P1a} \\ {{{subjectto}{{r_{j} = {\sum\limits_{i \in N}{\int_{0}^{1}{{\log\left( {1 + \frac{{p_{i\rightarrow j}(f)}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}{p_{l\rightarrow q}(f)}}}}}} \right)}f}}}},{\forall{j \in K}}}}} & {P1b} \\ {{\sum\limits_{j \in K}\left| {p_{i\rightarrow j}(f)} \middle| {}_{0}{\leq 1} \right.},{\forall{i \in N}},{\forall{f \in \left\lbrack {0,1} \right\rbrack}}} & {P1c} \\ {{0 \leq {p_{i\rightarrow j}(f)} \leq P_{\max}},{\forall{i \in N}},{j \in K},{\forall{f \in {\left\lbrack {0,1} \right\rbrack.}}}} & {P1d} \end{matrix}$

P1 can be difficult to solve directly due to two factors: (1) infinite problem size because of the continuous variable {p(f)}, and (2) the nonconvex constraint P1c. As discussed in more detail below, P1 is reformulated to a problem with a finite number of variables, and an efficient algorithm is then proposed to overcome the aforementioned difficulties.

The following proposition is described herein as Definition 1: A power allocation {p(f)} is said to be piecewise constant if the interval [0,1] can be partitioned into a finite number of sub-intervals such that, for every i ∈N and j ∈K, p_(i>1)(f) is constant on each one of those sub-intervals. Additionally, the following proposition is referred to as Definition 2: A piecewise constant power allocation is said to be k-sparse if it suffices to partition [0,1] into k sub-intervals such that the PSDs are flat on all of those sub-intervals.

The feasible rate region R can be defined as follows:

R={r ∈R ^(k): (P1b), (P1c), and (P1d) hold for all i ∈N, j ∈K, f ∈[0, 1]}.   Eq. 10:

Theorem 1, which is proven below, is as follows: The optimal utility of P1 can be attained by a k-sparse piecewise constant power allocation. Theorem 1 implies that any optimal rate vector r can be achieved by dividing the entire spectrum into k constant-PSD intervals with bandwidths β¹, . . . , β^(k). To reformulate P1, L={1, . . . , k} and power allocation over the j-th interval can be represented by a vector p^(j)=(p_(i→j) ^(j), . . . , p_(k→n) ^(j)). Therefore, P1 can reformulated as P2:

$\begin{matrix} {\mspace{79mu} {2\underset{p,\beta}{maximize}\mspace{14mu} {U(r)}}} & {P2a} \\ {{{{subjectto}\mspace{14mu} r_{j}} = {\sum\limits_{m \in L}{\beta^{m}{\sum\limits_{i \in N}{\log\left( {1 + \frac{p_{i\rightarrow j}^{m}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}^{m}}}}}} \right)}}}}},\mspace{79mu} {\forall{j \in K}}} & {P\; 2b} \\ \begin{matrix} {\mspace{79mu} {{\sum\limits_{j \in K}\left| p_{i\rightarrow j}^{m} \middle| {}_{0}{\leq 1} \right.},}} & {{\forall{i \in N}},{\forall{m \in L}}} \end{matrix} & {P2c} \\ {\mspace{79mu} {{{\sum\limits_{m \in L}\beta^{m}} = 1},}} & {P2d} \\ {\mspace{79mu} {{0 \leq p_{i\rightarrow j}^{m} \leq P_{\max}},{\forall{i \in N}},{j \in K},{\forall{m \in L}}}} & {P2e} \\ {\mspace{79mu} {{\beta^{m} \geq 0},{\forall{m \in {L.}}}}} & {P2f} \end{matrix}$

In P2, there are totally k+k²n=O(k²) variables. P2 is still difficult to solve because of the nonconvex constraints P2c. Moreover, for a large network that includes thousands of APs and UEs, P2 still has millions of variables, which prohibits the space and time complexities of solving P2. Described below are important facts and propositions that allow one to solve a simple case of P2 in which the utility function U(r) is an affine function of r. This technique can then be used to develop a scalable solution.

If the utility function U(r) is an affine function of r, one can arrive at the following proposition. Proposition 1: Suppose U(r) is an affine function of r. Then the maximum utility in P2 can be attained by letting each AP apply a flat power spectral density (PSD) over the entire spectrum (i.e., 1-sparse) to serve a single UE. To prove this proposition, the affine utility function can be written as:

$\begin{matrix} {{{U(r)} = {d + {\sum\limits_{j \in K}{c_{j}r_{j}}}}},} & {{Eq}.\mspace{14mu} 11} \end{matrix}$

for some constants d and (c_(j)). Therefore, P2 can be reformulated as P3, as follows:

$\begin{matrix} {{{2\underset{p,\beta}{maximize}\mspace{20mu} d} + {\sum\limits_{m \in L}{\beta^{m}{\sum\limits_{j \in K}{c_{j}{\sum\limits_{i \in N}{\log\left( {1 + \frac{p_{i\rightarrow j}^{m}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}^{m}}}}}} \right)}}}}}}},\mspace{20mu} {{subject}\mspace{14mu} {to}\mspace{14mu} \left( {P\; 2c} \right)},\left( {P\; 2d} \right),\left( {P\; 2e} \right),{{and}\mspace{14mu} {\left( {P\; 2f} \right).}}} & {P3a} \end{matrix}$

One can let (p, β) be a feasible solution to P3, and can define m* ∈L as a maximizer of

$\sum\limits_{j \in K}{c_{j}{\sum\limits_{i \in N}{{\log\left( {1 + \frac{p_{i\rightarrow j}^{m}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}^{m}}}}}} \right)}.}}}$

It can be seen that for any fixed power allocation p, the utility is always maximized by letting β^(m)*=1. That is, by letting each AP apply a flat power spectral density over the entire spectrum to serve only one UE.

Now, P3 can be reformulated to remove the l₀ norm constraint (P2c) and to further reduce the number of variables. The variable Φ_(i) can denote the set of users who can potentially be served by AP i, and the variable s_(i) can denote the UE to be served by AP i. Based on these variables and on Proposition 1 from above, P3 can be reformulated as P4 as follows:

$\begin{matrix} {2\underset{s,p}{maximize}{\sum\limits_{i \in N}{c_{s_{i}}{\log\left( {1 + \frac{p_{i}g_{i\rightarrow s_{i}}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{p_{l}g_{l\rightarrow s_{i}}}}}} \right)}}}} & {P4a} \\ {{{{subjectto}\mspace{14mu} 0} \leq p_{i} \leq P_{\max}},{\forall{i \in N}}} & {P4b} \\ {{s_{i} \in \Phi_{i}},{\forall{i \in {N.}}}} & {P4c} \end{matrix}$

In P4, s refers to the UE decisions for APs. Specifically, s_(i) is the index of the UE served by AP i, and it is possible for multiple APs to serve the same UE, i.e., s_(i)=s_(i′) for i≠i′. The variable s_(i)=Ø if AP i serves no UE, and one can define c_(Ø)=0. p=(p₁, . . . , p_(n)) as the vector that includes the constant PSDs of all APs over the entire spectrum. P4 has a total of 2n variables.

P4 can be viewed as a fractional programming problem, which has been well studied in prior works. A similar technique is used herein. First, it can readily be proven that P4 is equivalent to P5:

$\begin{matrix} {{2\underset{s,p,\gamma}{maximize}{\sum\limits_{i \in N}{c_{s_{i}}{\log \left( {1 + \gamma_{i}} \right)}}}} - {\sum\limits_{i \in N}{c_{s_{i}}\gamma_{i}}} + {\sum\limits_{i \in N}\frac{{c_{s_{i}}\left( {1 + \gamma_{i}} \right)}p_{i}g_{i\rightarrow s_{i}}}{n_{0} + {\sum\limits_{l \in N}{p_{l}g_{l\rightarrow s_{i}}}}}}} & {P5a} \\ {\mspace{79mu} {{{{subjectto}\mspace{14mu} 0} \leq p_{i} \leq P_{\max}},{\forall{i \in N}}}} & {P5b} \\ {\mspace{79mu} {{s_{i} \in \Phi_{i}},{\forall{i \in {N.}}}}} & {P5c} \end{matrix}$

When s and p are fixed, one can obtain the optimal γ as follows:

$\begin{matrix} {\gamma_{i}^{*} = {\frac{p_{i}g_{i\rightarrow s_{i}}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{p_{l}g_{l\rightarrow s_{i}}}}}.}} & {{Eq}.\mspace{14mu} 12} \end{matrix}$

The equivalence can be established by noting that the objective (P5a) is concave in (γ_(i))_(i ∈N) and that (P4a) is equal to (P5a) with γ in it. Then, based on a proven theorem, it can be shown that P5 is equivalent to P6:

$\begin{matrix} {{2\underset{s,p,\gamma,y}{maximize}\mspace{14mu} {\sum\limits_{i \in N}{c_{s_{i}}{\log \left( {1 + \gamma_{i}} \right)}}}} - {\sum\limits_{i \in N}{c_{s_{i}}\gamma_{i}}} + {\sum\limits_{i \in N}\left( {{2y_{i}\sqrt{{c_{s_{i}}\left( {1 + \gamma_{i}} \right)}p_{i}g_{i\rightarrow s_{i}}}} - {y_{i}^{2}\left( {n_{0} + {\sum\limits_{l \in N}{p_{l}g_{l\rightarrow s_{i}}}}} \right)}} \right)}} & {P6a} \\ {\mspace{79mu} {{{{subjectto}\mspace{14mu} 0} \leq p_{i} \leq P_{\max}},{\forall{i \in N}}}} & {P6b} \\ {\mspace{79mu} {{s_{i} \in \Phi_{i}},{\forall{i \in {N.}}}}} & {P6c} \end{matrix}$

When all other variables are fixed, one can obtain the optimal y:

$\begin{matrix} {y_{i}^{*} = {\frac{\sqrt{{c_{s_{i}}\left( {1 + \gamma_{i}} \right)}p_{i}g_{i\rightarrow s_{i}}}}{n_{0} + {\sum\limits_{l \in N}{p_{l}g_{l\rightarrow s_{i}}}}}.}} & {{Eq}.\mspace{14mu} 13} \end{matrix}$

Similarly, the optimal power allocation for other fixed other variables is:

$\begin{matrix} {p_{i}^{*} = {\min {\left\{ {P_{\max},\frac{{c_{s_{i}}\left( {1 + \gamma_{i}} \right)}g_{i\rightarrow s_{i}}y_{i}^{2}}{\left( {\sum\limits_{l \in N}{\gamma_{l}^{2}g_{i\rightarrow s_{l}}}} \right)^{2}}} \right\}.}}} & {{Eq}.\mspace{14mu} 14} \end{matrix}$

Finally, since the the utility function (P6a) can be decoupled in terms of each AP, the optimal UE decision variables s can be obtained through a simple search to find the user that maximizes the utility for each AP. This relation is depicted in Eq. 15 below:

$\begin{matrix} {s_{i} = \left\{ \begin{matrix} {Ø,{{{if}\mspace{14mu} {\max\limits_{j \in \Phi_{i}}\left\{ {{c_{j}\mspace{14mu} {\log \left( {1 + \gamma_{i}} \right)}} - {c_{j}\gamma_{i}} + {2y_{i}\sqrt{{c_{j}\left( {1 + \gamma_{i}} \right)}p_{i}g_{i\rightarrow j}}} - {y_{i}^{2}{\sum\limits_{l \in N}{p_{l}g_{l\rightarrow j}}}}} \right\}}} < 0}} \\ {{\arg \mspace{11mu} {\max\limits_{j \in \Phi_{i}}\left\{ {{c_{j}\mspace{11mu} {\log \left( {1 + \gamma_{i}} \right)}} - {c_{j}\gamma_{i}} + {2y_{i}\sqrt{{c_{j}\left( {1 + \gamma_{i}} \right)}p_{i}g_{i\rightarrow j}}} - {y_{i}^{2}{\sum\limits_{l \in N}{p_{l}g_{l\rightarrow j}}}}} \right\}}},{otherwise}} \end{matrix} \right.} & {{Eq}.\mspace{14mu} 15} \end{matrix}$

The proposed joint power control and user association strategy for the affine utility function is summarized in the algorithm below. It can be readily verified that this algorithm is guaranteed to converge, with the utility (P4a) monotonically non-decreasing after each iteration. The total number of variables needed to be solved at each iteration is O(n), moreover, and each iteration in the algorithm is performed in closed form. As discussed in more detail below, all of these propositions can be combined to solve P1 for any utility functions.

Algorithm 1: Joint Power Control and User Association for the Affine Utility Function

Initialization: Initialize s, p, γ.

Repeat:

-   -   Step 1: Update γ by Eq. 12;     -   Step 2: Update γ by Eq. 13;     -   Step 3. Update s by Eq. 15;     -   Step 4. Update p by Eq. 14;

Until Convergence.

As noted above, the optimal solution to P1 can be obtained by partitioning the entire spectrum into at most k segments, and by having each AP apply a flat PSD to serve only one UE on each segment. Motivated by this property, Algorithm 2 is designed to identify one segment with power allocation and user association in each iteration by optimizing the first order approximation of the utility function using Algorithm 1. The identified segment is added to the set of collected segments which is used to re-optimize P2 in each iteration. Since the set of collected segments grows after each iteration, the utility function increases after each iteration until it converges.

Algorithm 2: Interative Algorithm for Segment Pursuit

Initialization: Pick full spectrum reuse power allocation denoted by p⁰ as the initial feasible solution which gives a feasible rate vector of r⁰; t←0; P←{p⁰}.

Repeat:

-   -   Step 1: Solve P2 with the objective function being <∇U(r^((t))),         r> using Algorithm 1, which results in p^(t+1);     -   Step 2. If p^(t+1)∉P, P←P∪{p^(t+1)}, otherwise add an arbitrary         new segment that is not in P, t←t+1.     -   Step 3. Solve P2 given power allocations in P and obtain the         rate vector r^((t)).

Until Convergence.

Algorithm 2 can be interpreted as a Frank-Wolfe type Algorithm (also referred to as a conditional gradient algorithm). The main difference between Algorithm 2 and the conventional algorithm is that a line search is not performed. Rather, in each iteration of Algorithm 2, one segment in which each AP applies a flat PSD to serve at most one UE is found and added to the set P to re-optimize P2. The set P grows after each iteration because either a new segment is found or a random new segment is added. Based on Theorem 1, it is possible to find the global optimal solution for P1 after running Algorithm 2 for k iterations.

The performance gain of the proposed allocation scheme was investigated by comparing it with the following baseline schemes: (1) Full-spectrum-reuse +maximum reference signal receive power (MaxRSRP), in which every AP reuses all available spectrum with full transmit power and every UE is associated to the strongest AP in terms of the received power; (2) Full-spectrum-reuse+optimal user association, in which every AP reuses all available spectrum with full transmit power and user association is optimized for the utility; and (3) A pattern-pursuit algorithm which optimizes both spectrum allocation and user association assuming full transmit power at each AP. FIG. 2 depicts a table with parameter configurations used to test the performance gain of the proposed allocation scheme in accordance with an illustrative embodiment. The parameters in FIG. 2 are compliant with the LTE standard.

The performance of the proposed scheme was first tested in a network of medium size. FIG. 3 depicts a test deployment in a medium scale network in accordance with an illustrative embodiment. Specifically, 100 APs and 250 UEs were dropped over a 1330×1330 m² area as shown in FIG. 3. FIG. 4 depicts the objective of average packet delay versus traffic intensity curves for the deployment of FIG. 3 in accordance with an illustrative embodiment. From FIG. 4, it can be observed that as the average UE traffic increases to above 8, 13, and 30 packets/second, respectively, all three of the above-referenced baseline schemes fail to support all of the UEs. Conversely, the proposed solution has a significantly larger throughput (above 39 packets/second) as compared to the baseline schemes. The proposed solution also significantly reduces the delay, especially in the high traffic regime. The reason is that the proposed solution jointly optimizes the power allocation, spectrum allocation, and user association.

FIG. 5 depicts a test deployment in a large scale network in accordance with an illustrative embodiment. In this test scenario, the proposed allocation scheme is used to compute the allocation solution for a large network with the same density as the medium-scale network discussed above. The large network deployment depicted in FIG. 5 includes1000 APs and 2500 UEs over a 4200×4200 m² area. FIG. 6 depicts the objective of average packet delay versus traffic intensity curves for the deployment of FIG. 5 in accordance with an illustrative embodiment. As depicted in FIG. 6, the proposed solution has significantly larger throughput (above 30 packets/second) than full-spectrum reuse with maxRSRP association (7 packets/second), full-spectrum reuse with optimal user association allocation (14 packets/second), and centralized optimization with full power (22 packets/second). The proposed solution also outperforms baseline schemes in delay, especially in the high traffic regime.

Thus, described herein is an improved solution to the joint resource allocation problem in large networks. A scalable reformulation has been developed by exploiting the hidden sparse structure of the network, and an efficient algorithm is proposed to solve the reformulated problem with guranteed convergence. Moreover, each iteration is performed in closed form, which makes centralized resource allocation practically feasible even for a very large network. Numerical results show substantial gains compared to all the other baseline schemes for networks with 1000 access points or more. The proposed algorithm applies to any concave utility functions such as sum rate, minimum user service rate (max-min fairness), and sum log-rate (proportional fairness).

FIG. 7 is a block diagram depicting a central controller 700 for the proposed system in accordance with an illustrative embodiment. The central controller 700 includes a processor 705, a memory 710, a transceiver 715, and an interface 720. In alternative embodiments, the central controller 700 may include additional, fewer, and/or different components. The memory 710 can be any type of computer-readable medium that is configured to store computer-readable instructions to implement any of the operations, programs, and algorithms described herein. The processor 705 is configured to execute the computer-readable instructions stored in the memory 710. The processor 705 is also configured to control the transceiver 715 such that packets and other information can be received and transmitted by the central controller 700. The transceiver 715, which is used to receive/transmit information, can be any type of transceiver. The transceiver 715 can be in the form of a single component that both sends and receives, or alternatively the transceiver 715 can be in the form of a distinct receiver and a distinct transmitter. The interface 720 can be any type of input/output (I/O) component such as a display, keyboard, mouse, data port, etc. that allows an individual or another computing device to interact with the central controller 700.

Proof of Theorem 1

The above-discussed Theorem 1 can be proven by first proving that any point in the feasible rate region R given by Eq. 10 can be obtained with a (kn+1)-sparse piecewise constant power allocation. It is then further proven that a k-sparse piecewise constant power allocation suffices. For a network consisting of k UEs and n APs, there are totally kn links. The instantaneous spectral efficiency set S_(link) and the achievable rate region of all links R_(link) can be defined as follows:

$\begin{matrix} {S_{link} = {\left\{ {{{s \in {R^{kn}:s_{{{\langle{i - 1})}k} + j}}} = {\log\left( {1 + \frac{p_{i\rightarrow j}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}}}}}} \right)}},{\sum\limits_{j \in K}\left| p_{i\rightarrow j} \middle| {}_{0}{\leq 1} \right.},{0 \leq p_{i\rightarrow j} \leq P_{\max}},{\forall{i \in N}},{j \in K}} \right\} \mspace{14mu} {and}}} & {{Eq}.\mspace{14mu} 16} \\ {{R_{link} = \left\{ {{{r_{{kn} \times 1}:r_{i\rightarrow j}} = {\int_{0}^{1}{{\log\left( {1 + \frac{{p_{i\rightarrow j}(f)}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}{p_{l\rightarrow q}(f)}}}}}} \right)}f}}},{\sum\limits_{j \in K}\left| {p_{i\rightarrow j}(f)} \middle| {}_{0}{\leq 1} \right.},{0 \leq {p_{i\rightarrow j}(f)} \leq P_{\max}},{\forall{i \in N}},{j \in K},{f \in \left\lbrack {0,1} \right\rbrack}} \right\}}} & {{Eq}.\mspace{14mu} 17} \end{matrix}$

The instantaneous spectral efficiency set S_(link) is the set of all spectral efficiencies achievable by power allocation at any frequency. By the standard convexity argument, it is apparent that the achievable rate region of all links R_(link) is the convex hull of the instantaneous spectral efficiency set S_(link) i.e., R_(link)=conv{S_(link)}. Then by the Caratheodory theorem, it can be proven that the number of elements in S_(link) that is needed to construct an arbitrary data rate vector r_(kn×1)=(r_(i→j), ∀i ∈N, ∀j ∈K) in R_(link) is at most kn+1. In other words, any achievable rate vector can be obtained with at most kn+1 piecewise constant power allocations. That is, for any (r_(i→j))_(i ∈N, j ∈K) ∈R_(link), it follows that:

$\begin{matrix} {{r_{i\rightarrow j} = {\sum\limits_{m = 1}^{{kn} + 1}{\beta^{m}{\log\left( {1 + \frac{p_{i\rightarrow j}^{m}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}^{m}}}}}} \right)}}}},} & {{Eq}.\mspace{14mu} 18} \end{matrix}$

for some piecewise constant power allocations {p_(i→j) ^(m)}_(i ∈N, j ∈K, m=1,2, . . . , kn+1) and corresponding bandwidths {β^(m)}m=_(1,2, . . . , kn+1) with Σ_(m−1) ^(kn+1)β^(m)=1.

Next, the number of piecewise constant power allocations are further reduced to obtain an arbitrary achievable rate vector. Let r_(k/1)=(r_(j), ∀j ∈K) be an element in the set R. By definition, there exists a vector r′_(kn×1) ∈R_(link) such that

${r_{j} = {\sum\limits_{i \in N}r_{i\rightarrow j}^{\prime}}},{\forall{j \in {K.}}}$

Based on this portion of the proof, it is known that the vector r′ can be constructed using kn+1 piecewise constant power allocations denoted as {p_(i→j) ^(m)}_(i∈N, j∈K, m=1,2, . . . , kn+1). Then for each j ∈K, it can be seen that:

$\begin{matrix} {\begin{matrix} {r_{j} = {\sum\limits_{i \in N}r_{i\rightarrow j}^{\prime}}} \\ {= {\sum\limits_{i \in N}{\sum\limits_{m = 1}^{{kh} + 1}{\beta^{m}\log}}}} \\ {\left( {1 + \frac{p_{i\rightarrow j}^{m}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}^{m}}}}}} \right)} \\ {= {\sum\limits_{m = 1}^{{kh} + 1}{\beta^{m}{\sum\limits_{j \in N}\log}}}} \\ {\left( {1 + \frac{p_{i\rightarrow j}^{m}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}^{m}}}}}} \right)} \end{matrix}{A\mspace{14mu} {new}\mspace{14mu} {set}\mspace{14mu} {can}\mspace{14mu} {be}\mspace{14mu} {defined}\mspace{14mu} {as}\mspace{14mu} {follows}\text{:}}} & \begin{matrix} {{Eq}.\mspace{14mu} 19} \\ \; \\ {{Eq}.\mspace{14mu} 20} \\ \; \\ \; \\ \; \\ \; \\ \; \\ {{Eq}.\mspace{14mu} 21} \\ \; \end{matrix} \\ {T = \left\{ {{{t_{k \times 1}:t_{j}} = {{\sum\limits_{i \in N}{{\log\left( {1 + \frac{p_{i\rightarrow j}g_{i\rightarrow j}}{n_{0} + {\sum\limits_{l \in {N:{l \neq i}}}{g_{l\rightarrow j}{\sum\limits_{q \in K}p_{l\rightarrow q}}}}}} \right)}{\sum\limits_{j \in K}{p_{i\rightarrow j}}_{0}}}} \leq 1}},{0 \leq p_{i\rightarrow j} \leq P_{\max}},{\forall{i \in N}},{j \in K}} \right\}} & {{Eq}.\mspace{14mu} 22} \end{matrix}$

It is apparent that the vector r_(k×1), is in the convex hull of T. Therefore, by the Caratheodory theorem, it can be proven that the number of elements in T that is needed to construct an arbitrary data rate vector r ∈R is at most k+1. Moreover, r is either on the boundary of the convex hull of T or it is an interior point. In either case, there exists r* on the boundary that dominates r in every dimension. Therefore, the number of elements in T that is needed to construct an optimal rate vector r_(k×1) in R is at most k, which implies that an optimal rate vector r_(k×1) can be achieved by dividing the spectrum [0,1] into k intervals with lengths β₁, . . . , β_(k) and assigning power vectors p¹, . . . , p^(k) to these intervals, where p^(m)={p_(i→j) ^(m)}_(i ∈N, j ∈K).

In an illustrative embodiment, any of the operations described herein may be performed by a computing system that includes a memory, processor, user interface, transceiver, and any other computing components. The operations can be stored as computer-readable instructions on a computer-readable medium such as the computer memory. Upon execution by the processor, the computer-readable instructions are executed as described herein.

The word “illustrative” is used herein to mean serving as an example, instance, or illustration. Any aspect or design described herein as “illustrative” is not necessarily to be construed as preferred or advantageous over other aspects or designs. Further, for the purposes of this disclosure and unless otherwise specified, “a” or “an” means “one or more”.

The foregoing description of illustrative embodiments of the invention has been presented for purposes of illustration and of description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the invention. The embodiments were chosen and described in order to explain the principles of the invention and as practical applications of the invention to enable one skilled in the art to utilize the invention in various embodiments and with various modifications as suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims appended hereto and their equivalents. 

What is claimed is:
 1. A system for allocating resources in a communications network, the system comprising: a plurality of access points, wherein each access point in the plurality of access points is in communication with one or more user devices; and a central controller in communication with the plurality of access points, wherein the central controller is configured to: receive traffic information and channel information for access point—user device links in the communications network; receive an intensity of independent homogeneous Poisson traffic intended for the one or more user devices; and determine a joint resource allocation based at least in part on the received traffic and channel information and the intensity of independent homogeneous Poisson traffic, wherein the joint resource allocation is determined with an algorithm.
 2. The system of claim 1, wherein the central controller is further configured to receive interference information for the access point—user device links in the communications network.
 3. The system of claim 2, wherein the central controller determines the joint resource allocation based at least in part on the interference information.
 4. The system of claim 1, wherein the central controller is further configured to transmit the joint resource allocation to the plurality of access points.
 5. The system of claim 4, wherein the joint resource allocation comprises spectrum allocation, user association, and power allocation.
 6. The system of claim 1, wherein the plurality of access points are configured to implement the joint resource allocation.
 7. The system of claim 6, wherein each of the plurality of access points is configured to apply a continuous power spectral density over a frequency band used in the communications network.
 8. The system of claim 1, wherein the algorithm is derived based on a non-convex infinite-size joint spectrum allocation, user association, and power allocation problem.
 9. The system of claim 8, wherein the algorithm is based on a reformulation of the non-convex infinite-size joint spectrum allocation, user association, and power allocation problem into an optimization problem with finite size.
 10. The system of claim 1, wherein the central controller is configured to model average channel conditions based on path loss and statistics of small scale fading.
 11. The system of claim 1, wherein the central controller is further configured to determine a power spectral density for each of the access point—user device links over an entire frequency spectrum used by the communications network.
 12. The system of claim 11, wherein the central controller determines the joint resource allocation based at least in part on the determined power spectral densities.
 13. A method for allocating resources in a communications network, the method comprising: receiving, by a central controller and from a plurality of access points in communication with the central controller, traffic information and channel information for access point—user device links in a communications network, wherein each access point in the plurality of access points is in communication with one or more user devices; receiving, by the central controller, an intensity of independent homogeneous Poisson traffic intended for each of the one or more user devices; and determining, using an algorithm, a joint resource allocation based at least in part on the received traffic and channel information and the intensity of independent homogeneous Poisson traffic.
 14. The method of claim 13, further comprising receiving interference information for the access point—user device links in the communications network.
 15. The method of claim 14, further comprising determining, by the central controller, the joint resource allocation based at least in part on the interference information.
 16. The method of claim 13, further comprising transmitting the joint resource allocation to the plurality of access points.
 17. The method of claim 16, wherein the joint resource allocation comprises spectrum allocation, user association, and power allocation.
 18. The method of claim 13, further comprising implementing, by the plurality of access points, the joint resource allocation.
 19. The method of claim 18, wherein implementing the joint resource allocation comprises, in part, applying a continuous power spectral density over a frequency band used in the communications network.
 20. The method of claim 13, further comprising determining, by the central controller, a power spectral density for each of the access point—user device links over an entire frequency spectrum used by the communications network. 